clear all;
x0=[3,-1,0,1]';
eps = 1.0e-4;
x_star = Newton(x0,eps); %阻尼牛顿法
x_star = FD(x0,eps);     %最速下降
x_star = DFP(x0,eps);    %DFP方法
x_star = FR(x0,eps);     %FR法
